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1.1 Introduction 

This chapter is intended to give an introduction to the theory of thermal 
fluctuations and their consequences for static and dynamic correlations in ge- 
ometrically frustrated antiferromagnets, focusing on the semiclassical limit, 
and to discuss how our theoretical understanding leads to an explanation of 
some of the main observed properties of these systems. A central theme will 
be the fact that simple, classical models for highly frustrated magnets have a 
ground state degeneracy which is macroscopic, though accidental rather than 
a consequence of symmetries. We will be concerned in particular with: (i) 
the origin of this degeneracy and the possibility that it is lifted by thermal 
or quantum fluctuations; (ii) correlations within ground states; and (iii) low- 
temperature dynamics. We concentrate on Heisenberg models with large spin 
S, referring to the chapter by G. Misguich for a discussion of quantum spin 
liquids, and to the chapter by M. Gingras for an overview of geometrically 
frustrated Ising models in the context of spin ice materials. Several earlier re- 
views provide useful further reading, including [1], [2] and [3] for experimental 
background, and [4] and [5] for an alternative perspective on theory. 

To provide a comparison, it is useful to begin by recalling the behaviour 
of an unfrustrated antiferromagnet. To be definite, consider the Heisenberg 
model with nearest neighbour exchange J on a simple cubic lattice. As the 
lattice is bipartite - it can be separated into two interpenetrating sublattices, 
in such a way that sites of one sublattice have as their nearest neighbours 
only sites from the other sublattice - the classical ground states are two- 
sublattice Neel states, in which spins on one sublattice all have the same 
orientation, and those on the other sublattice have the opposite orientation. 
These states are unique up to global spin rotations, which are a symmetry of 
the model. Their only low energy excitations are long wavelength spinwaves. 
These are Goldstone modes - a consequence of the symmetry breaking in 
ground states - and have a frequency u>(k) that is linear in wavevector k at 
small k. This classical picture carries over to the quantum system, and for 
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S 3> 1 it is sufficient to take account of fluctuations using harmonic spinwave 
theory. In particular, within this approximation the sublattice magnetisation 
at low temperature is reduced from its classical ground state value S by an 
amount 



where {n(k)) is a Bose factor giving the number of thermally excited spin 
waves at wavevector k and 1/2 represents the zero-point contribution, with 
the integral running over the Brillouin zone of volume fl. Fluctuations increase 
with temperature and the sublattice magnetisation falls to zero at the Neel 
temperature Tn- Within mean field theory IcbTn = zJS 2 , where z is the 
number of nearest neighbour sites (six on the simple cubic lattice). 

A central reason for the interest in geometrically frustrated magnets is that 
they hold out the possibility of evading Neel order. At the simplest level, there 
is a tendency for frustrated systems to have many low-frequency modes, which 
means both that excitations are effective in reducing the ordered moment, 
because of the factor of zJS/fiu>(\£) in and that they are thermally 

populated even at low temperature. More fundamentally, we will see that 
frustration may lead to classical ground states of a quite different kind and 
suppress Tn to zero. 

Since the term frustration is used in several different contexts, it is worth- 
while to set out some distinctions before going further. In general terms, 
classical frustrated systems have Hamiltonians with competing interactions 
which make contributions to the energy that cannot simultaneously be min- 
imised. The concept was originally discussed in relation to spin glasses, but 
these are set apart from the systems we are considering here by the fact that 
quenched disorder has a controlling influence on their properties. Frustration 
as a way of destabilising Neel order has been studied extensively in models 
with competing nearest neighbour and further neighbour interactions, notably 
the Ji — J2 model on the square lattice [6], illustrated in Fig. 11.11 The classi- 
cal ground state of this model depends on the ratio J1/J2. For Ji > 2J 2 > 
neighbouring spins are antiparallel, enforcing ferromagnetic alignment of sec- 
ond neighbours and frustration of the interaction J 2 . In the other regime, 
2.J2 > J\ > 0, second neighbours are antiferromagnetically aligned at the 
expense of frustration of half of the Ji interactions. Interest focuses on the 
point Ji = 2J2 where these alternative classical states are degenerate, and the 
consequences of frustration are likely to be largest. While models of this kind 
provide an attractive starting point for theoretical work, there are likely to 
be difficulties in finding experimental realisations with interaction strengths 
that place them close to the degeneracy point. From this perspective, the 
long-appreciated [7, 8] attraction of geometrically frustrated magnets is that 
they are systems in which structure alone may destabilise Neel order, with 
only nearest neighbour interactions. To illustrate this at an elementary level, 
consider Heisenberg spins at the vertices of two corner-sharing triangles with 
nearest-neighbour antiferromagnetic interactions, also shown in Fig. 11.11 The 
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Fig. 1.1. Left and centre: Ji — J2 model, showing ground state spin configurations 
for: 2Ji > J2 (left); and J2 > 2Ji (centre). Right: ground states of classical Heisen- 
berg spins at vertices of two corner-sharing triangles, with degeneracy arising from 
rotations about the common spin, as indicated. 

ground states are configurations in which spins within each triangle are copla- 
nar and at relative angles of 2ir/3. They have an accidental degeneracy (in 
addition to that arising from symmetry) under relative rotations of the spin 
planes for the two triangles about the axis defined by the orientation of their 
common spin. 

1.2 Models 

The models we are concerned with extend some features present in the sim- 
ple system of two corner-sharing triangles to a periodic lattice. In general we 
will consider non-bipartite lattices constructed from corner-sharing arrange- 
ments of frustrated clusters, with local magnetic moments at the vertices of 
each cluster and exchange interactions of equal strength between all moments 
in each cluster (other arrangements are also of interest, but typically show 
less dramatic consequences of frustration) [9, 10]. An important example in 
two dimensions is the kagome lattice, formed from corner-sharing triangles; a 
three-dimensional analogue is the pyrochlore lattice, built from corner-sharing 
tetrahedra: see Fig. 11.21 for illustrations of both. 




Fig. 1.2. Left: kagome lattice. Right: pyrochlore lattice 



The Hamiltonian for these models, written in terms of the exchange energy 
J and the spin operators at sites i, has the form 
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H = J^2s i 'S j = -^2\L a \ 2 + c 

{ij) a 

where L a = ^ S, . (1.2) 

Here the first sum runs over neighbouring pairs ij of sites, while the second 
sum is over clusters a. To recognise that this second expression is a correct 
rewriting of 7i in terms of the total spin L a of each cluster a, it is necessary 
only to note that expansion of |L Q | 2 generates the required off-diagonal terms 
Si • Sj, together with diagonal terms S 2 that contribute to the constant c. 
The apparent simplicity of this second form is of course deceptive, since the 
operators L Q and associated with two clusters a and (3 that share a site 
are not independent. 

For future reference it is useful to introduce some terminology. The frus- 
trated clusters are in general simplices, and their centres occupy sites of a sec- 
ond lattice, called the simplex lattice. Spins in our models arc located at the 
mid-points of nearest-neighbour links of the simplex lattice. For the kagome 
magnet the simplex lattice is the honeycomb lattice, and for the pyrochlorc 
magnet it is the diamond lattice. 

While the Hamiltonian of (|1.2[) provides a useful basis for understand- 
ing the properties of a range of geometrically frustrated magnetic materi- 
als, various additional physical contributions to a realistic model may also 
be important. These include single- ion anisotropy [12], further neighbour ex- 
change [13, 14], dipolar interactions [15, 16], Dzyaloshinskii-Moriya interac- 
tions [17], magnetoelastic coupling [18], site dilution [19,20] and exchange 
randomness [21]. In many cases the associated energy scales are small. They 
set a temperature scale much smaller than nearest-neighbour exchange, below 
which they may induce magnetic order or spin freezing, but they can be ne- 
glected at higher temperatures. We omit all these perturbations and restrict 
our discussion to models with only nearest neighbour exchange. 

1.3 Some Experimental Facts 

The single most revealing property of a geometrically frustrated magnet is 
arguably the dependence on temperature T of its magnetic susceptibility \- 
It is convenient to consider plots of x^ 1 vs T, which at high temperature have 
the linear form 

X - x c<T-0 cw , (1.3) 

where the Curie- Weiss constant Oqw characterises the sign and strength of 
interactions. In an antiferromagnet <9cw is negative, and for the model of 
(|1.2[) one has ke6>cw = —zJS 2 . Without frustration, magnetic order, sig- 
nalled by a cusp in Xi appears below the Neel temperature, Tn ~ |6>cw|- 
By contrast, in geometrically frustrated systems nothing sharp is observed at 
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Fig. 1.3. Characteristic behaviour of a geometrically frustrated antiferromagnet. 
Left: sketch of x~ vs T. Right: sketch of S(Q,u>) vs Q 

the temperature scale set by interaction strength: instead, the paramagnetic 
phase extends to temperatures T <c 6>cw- Ordering or spin freezing may ap- 
pear at a lower temperature T c , but a large value for the ratio / = \0cw\/T c 
is a signature of frustration [1]. This behaviour is illustrated schematically in 
Fig. 11.31 references to experimental papers are given in Table [TTT1 

More detailed information on low temperature behaviour is provided by 
magnetic neutron scattering (see the chapter by S. T. Bramwell). Again, we 
sketch typical observations in Fig. 11.31 and give references in Table 11.11 The 
dynamical structure factor S(Q, to) has a broad peak at finite wavevector Q, 
showing that spin correlations are predominantly short-range and antiferro- 
magnetic. The width of this peak indicates a correlation length of order the 
lattice spacing, while the small value of the elastic scattering cross-section for 
Q — > shows that correlations suppress long wavelength fluctuations in mag- 
netisation density. This form stands in contrast both to that in unfrustrated 
antiferromagnets, where Neel order leads to magnetic Bragg peaks, and to 
that in systems with short-range ferromagnetic correlations, where the struc- 
ture factor is peaked at Q = 0. Inelastic scattering has a width in frequency 
u) that decreases with decreasing temperature, and in materials that show 
spin freezing, scattering weight is transferred from the inelastic to the elastic 
response with little change in Q-dependence on cooling through T c . 



Table 1.1. Three geometrically frustrated antiferromagnets 



Material 


Structure 


|G>cw| 


T c 


References 


SrGa 3 Cr 9 0i 9 


pyrochlore slabs 


515 K 


4 K 


[22-25] 


hydronium iron jarosite 


kagome 


700 K 


14 K 


[26] 


Y2M02O7 


pyrochlore 


200 K 


22 K 


[27, 28] 
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Properties of three well-studied geometrically frustrated antiferromagnets 
are set out in Table II. JFl Two basic theoretical questions arise. Why is there 
no magnetic ordering at T ~ |6>cw|? And what is the nature of correlations 
in the strongly interacting regime T <C 6>cw? 

1.4 Classical Ground State Degeneracy 

To get insight into the answers to these questions, we start by considering 
ground states of models defined by (|1.2[) in the classical limit, in which the 
Si are not operators but three-component vectors of magnitude S. As a first 
step, it is useful to examine a single tetrahedral cluster of four spins, with the 
Hamiltonian 

H = ^\L\ 2 + c where L = Si + S 2 + S 3 + S 4 . (1.4) 

By writing the Hamiltonian in terms of the cluster spin L we see at once that 
ground states are those with L = 0. Such an arrangement of four vectors, 
each having three components, with resultant zero is shown in Fig. 11.41 these 
ground states have two internal degrees of freedom, indicated in Fig. 11.41 by 
the angles 9 and <j>, in addition to the degeneracies under global rotations 
which are expected from the symmetry of TL. 




Fig. 1.4. A ground state configuration for a frustrated cluster of four classical 
Heisenberg spins 

We should next understand how this accidental ground state degeneracy 
extends from a single cluster to a periodic lattice. We can do so using a 
counting argument [9,10], which compares F, the number of degrees of freedom 
in the system with K, the number of constraints that must be satisfied in 
ground states. The central point is that if all constraints are independent, 
then the number of ground state degrees of freedom is given by the difference 



All three examples show spin freezing below T c 
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F — K. Such an argument was used by Maxwell in 1864 to discuss the stability 
of mechanical systems of jointed rods [11], and is sometimes referred to as a 
Maxwellian counting argument. For a system of N s classical Heisenberg spins, 
F = 2N S , since two angles are required to specify the orientation of each 
spin. And in a system with the Hamiltonian of (|1.2[) consisting of N c clusters, 
K = 3N C , since in ground states all three components of L Q must be zero for 
every cluster a. Under the assumptions that all constraints can be satisfied 
simultaneously, and that they are all linearly independent, we arrive at an 
estimate for D, the number of ground-state degrees of freedom: D = F — K . 
Taking the example of the pyrochlore lattice, we have N s = 2N C (since four 
spins are associated with each tetrahedron, but every spin is shared between 
two tetrahedra) and hence D — N c , an extensive quantity. 

This is a striking conclusion: it suggests that there are local degrees of free- 
dom which can fluctuate independently without the system leaving its ground 
state manifold. The argument has two implications for our understanding of 
the experimental results summarised in Sect. 11.31 First, macroscopic degen- 
eracy may prevent long range order at the temperature scale set by interac- 
tion strength, since there are many low-energy configurations that lack order. 
Second, since the magnetisation of each cluster is zero in all ground states, 
the amplitude of long wavelength fluctuations in the magnetisation density is 
small at low temperature, and so the dynamical structure factor S(Q,u>) is 
small at low Q. 




Fig. 1.5. Illustration of how ground state degrees of freedom arise for the Heisenberg 
model on the kagome lattice: spins on the central hexagon may be rotated together 
through any angle about the axis defined by the outer spins, without leaving the 
ground state. 

At this point it is worth pausing to consider possible limitations to the 
counting argument that has been presented. As noted, it rests on an assump- 
tion that all ground state constraints are linearly independent. If this is not 
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the case, we underestimate D. In our context, corrections are important if 
they make an extensive contribution to D. This occurs in the kagome lattice 
Heisenberg antiferromagnet: in this case our estimate yields D — (since, 
for a lattice built from corner-sharing triangles, N s = 3N c /2), but by explicit 
construction one finds sets of states with special spin arrangements [29,30] 
for which D = N s /9. Such an arrangement is illustrated in Fig. 11.51 By con- 
trast, for the pyrochlore Heisenberg antiferromagnet, it is known [9, 10] that 
corrections to the estimate for D are at most sub-extensive. 

The view of classical geometrically frustrated Heisenberg antiferromagnets 
that emerges at this stage is summarised by the cartoon of phase space given 
in Fig. 11.61 within (high-dimensional) phase space for the system as a whole, 
the ground states form a manifold with a dimension that is much smaller but 
nevertheless extensive. At temperatures small compared to the Curie- Weiss 
constant (keT <C JS 2 ), the system is confined to a region of phase space that 
forms a thin layer around the ground state manifold. Quantum effects can be 
neglected provided JS <C keT, and so a strongly correlated, classical window, 
JS < k B T < JS 2 , opens for large S. 




1.5 Order by Disorder 

The fact that extensive ground state degeneracy in classical, geometrically 
frustrated antiferromagnets is, in the technical sense, accidental prompts us 
to ask whether it has robust consequences in the presence of thermal or quan- 
tum fluctuations. Specifically, since the degeneracy is not a consequence of 
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symmetry, one expects the spectrum of fluctuations around each ground state 
to be different: the possibility arises that ground states with the lowest ex- 
citation frequencies are selected, because they have the largest entropy and 
the smallest zero-point energy. Such an apparently paradoxical mechanism, by 
which fluctuations enhance order instead of suppressing it, is termed 'order- 
by-disorder' [31,32]. 

We will consider first the effects of thermal fluctuations, and begin by 
discussing a cluster of four spins. Two ground states with fluctuations of 
contrasting types are illustrated in Fig. 11.71 For the configuration shown on 
the left, the total spin of the cluster has a magnitude |L| that varies with the 
departure 59 from the ground state as |L| cx 59. Since the excitation energy 
is proportional to |L| 2 , it has a conventional, quadratic dependence on 59. By 
contrast, for the excitation from a collinear ground state shown on the right, 
|L| cx (59) 2 : this mode is therefore soft, with an energy proportional to (59) 4 . 
We wish to understand whether the presence of this soft mode leads almost 
collinear configurations to dominate at low temperature. 




Fig. 1.7. Fluctuations away from ground state configurations for a cluster of four 
spins. Left: a conventional fluctuation; right: a soft mode. 



Analysis for a cluster of four spins is simple enough that it can be followed 
through in full. To illustrate the range of possible outcomes, we will consider 
spins with n components, comparing behaviour for n = 3 and n = 2. We use 
the coordinate system shown in Fig. 11.81 Our aim is to evaluate the thermal 
probability distribution P n (a) of the angle a between the pair of spins Si and 
S2- The distribution P n (a)da is a product of two factors. One stems from the 
measure for S2, and is sin(a)da or da, for n — 3 or n — 2 respectively. The 
other comes from integrating over orientations of S3 and S4: it is 

Z n (a) cx JdS 3 J dS 4 exp (-^ |S 3 + S 4 - 2Scos(a/2)z| 2 ^ . (1.5) 

In the low temperature limit, this can be evaluated by expanding the en- 
ergy to quadratic order in deviations from a ground state. For Heisenberg 
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spins (n — 3) the low energy configurations have |— y | , \S\ <C 1; for n = 2, 
spins are coplanar and two coordinates are fixed: <p = 5 — 0. In a quadratic 
approximation the energy is 

— {(/3-7) 2 cos 2 (a/2) +[(/? + 7 ) 2 + £ 2 ] sin 2 (a/2)} , (1.6) 

so that (including for n = 3 a factor of sin 2 (a/2) arising from dSadSzt) 

Z 3 (a) cx [cos(a/2)] _1 and Z 2 cx [cos(a/2) sin(a/2)]~ 1 . (1.7) 
Combining contributions, we have 

P 3 (a) cx sin(a/2) and P 2 (a) cx — ^— . (1.8) 

sm(a) 

In this way we discover contrasting behaviour for the two cases. With n = 
3, the system explores all values of a even in the low temperature limit. But 
for n = 2 our unnormalised result for Zi(a) has non-integrable divergences 
at a — and a = 7r: in a more detailed treatment, retaining contributions to 
the energy quartic in coordinates, these divergences are cut off on a scale set 
by temperature, but in the low-temperature limit P2(a) approaches a sum of 
two delta functions, located at a = and a = it. Thus order by disorder is 
absent for n = 3 but perfect for n — 2. 




Passing from a single cluster to an extended system, consider the sketch 
of phase space given in Fig. 11.91 Here, repeating the convention of Fig. 11.61 
the shading indicates the region accessible at low temperature. One part of 



1 Fluctuations in geometrically frustrated magnets 



11 



this region is concentrated near points on the ground state manifold at which 
there are soft modes, as represented by in Fig. 11.91 by a bulge, while another 
part is distributed in the neighbourhood of the remainder of the ground state 
manifold. To decide whether the system displays order by disorder, we need 
to understand which of these two parts dominates. Introducing coordinates x 
and y, respectively parallel and perpendicular to the ground state manifold, 
in the low temperature limit we obtain a measure -P(x) on the ground state 
manifold by integrating over transverse fluctuations [9, 10]. Characterising 
these fluctuations by dynamical frequencies oji (x) , we obtain 



The system has extra soft modes (in addition to those associated with the 
ground state coordinates x) at points Xo where one or more of the harmonic 
frequencies wj(xo) vanishes. At these points -P(xo) is divergent. As for a cluster 
of four spins, behaviour depends on whether any such divergence is integrable. 
If it is, the system explores the whole of the ground state manifold in the low 
temperature limit, but if it is not, then those ground states with soft modes are 
selected by thermal fluctuations. Detailed considerations, tested using Monte 
Carlo simulations, show for the Heisenberg antiferromagnet on the kagome 
lattice that there is coplanar spin order in the low temperature limit [29], 
while on the pyrochlore lattice there is no order by disorder [9, 10]. 



Fig. 1.9. Schematic view of phase space. The full curve represents the ground state 
manifold. Coordinates x and y are respectively parallel and perpendicular to it. 

The possibility of ground state selection due to quantum fluctuations can 
be discussed using an approach similar in spirit to the one we have taken for 
thermal fluctuations, although the outcome has significant differences. Refer- 




(1.9) 
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ring again to Fig. 11.91 one can treat excitations around a particular point x on 
the ground state manifold using harmonic spin wave theory. Excitations in- 
volving the coordinates y locally orthogonal to the ground state manifold are 
conventional modes with non-zero frequencies u>i (x) , which have already made 
an appearance in (|1.9p . By contrast, fluctuations involving the coordinates x 
are, within a harmonic approximation, zero modes. The zero-point energy 
of the conventional, finite-frequency modes provides an effective Hamiltonian 
for these remaining degrees of freedom, the classical ground state coordinates. 
This Hamiltonian takes the form 



The components of x consist of pairs that are, within the approximations of 
harmonic spin wave theory (see (11.14| ). canonically conjugate. Treating them 
as classical commuting variables, the ground state is the set of points xq on 
which 7i c ff(x) is minimised. More accurately, the ground state wavefunction 
for large S is peaked at xg, but has zero-point fluctuations in an effective 
potential defined by W e ff( x )- 

It is not straightforward to anticipate what features of a classical ground 
state spin configuration will minimise 7i c g(x): since all w;(x) contribute, one 
could equally imagine focusing on either the highest frequencies or the lowest 
ones. In the examples that have been studied in detail, however, it seems that 
minima lie at points xg where some w/(x) vanish, which we have seen are 
also the states favoured by thermal fluctuations. In particular, for Heisenberg 
antiferromagnets at large S the selected spin configurations are coplanar on 
the kagome lattice [33] and collincar on the pyrochlore lattice [34]. In both 
examples, one third of the w;(x) become soft modes at the corresponding 
points xg: the coplanar or collinear configurations, respectively. 

The principal difference between the ordering effects of thermal and quan- 
tum fluctuations is that in the first case, as we have seen, order may or may 
not arise on the limit J>T, depending on the nature of the thermal ground 
state distribution -P(x), while in the second case we always expect order for 
S ^> 1, because by taking S sufficiently large, one can ensure that quantum 
fluctuations around the minimum of 7Y c ff (x) are arbitrarily small. Within this 
framework, the scenario by which one arrives at a spin liquid on reducing 
S is clear, at least in principle. For smaller S, the quantum fluctuations are 
larger and the ground state wavefunction is less well localised around the 
minimum of H c s (x), while below a critical value of S, the quantum ground 
state wavefunction becomes delocalised over the entire classical ground state 
manifold and the system loses magnetic order. At large, fixed S long range or- 
der induced by quantum fluctuations is suppressed thermally above a critical 
temperature T c ~ JS. While the expectation that spin liquids are favoured at 
small S is common to our discussions of both (1.1) and (1.10), one should of 
course remember that the two equations embody different physics: harmonic 
and anharmonic fluctuations, respectively. 




(1.10) 
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Efforts to identify experimental examples of order by disorder must face the 
problem of establishing that fluctuations, rather than additional interaction 
terms in the Hamiltonian, are the cause of what is observed. For the garnet 
Ca3Fe2Ge30i2, a material with two interpenetrating magnetic lattices coupled 
via zero-point fluctuations, it has been shown that a spinwave gap in the Neel 
ordered state indeed arises mainly in this way, by independent determination 
of the size of single ion anisotropy (the other possible origin for the gap) [35], 
and via the characteristic temperature dependence of the gap [36] . 

1.6 Ground State Correlations 

As we have seen, in some circumstance a model geometrically frustrated mag- 
net (for example, the classical Heisenberg model on the pyrochlore lattice) 
explores its entire ground state manifold at low temperature, and the exper- 
imental evidence from elastic and inelastic neutron scattering suggests that 
this is a reasonable picture for the behaviour of a range of frustrated magnetic 
materials. We are led to ask in this section whether there are any important 
correlations within the ground state manifold. We will find (for a large class of 
models: those in which the simplex lattice is bipartite) that there are indeed 
long-range correlations within ground states [10], and that these can be char- 
acterised in terms of fluctuations of a Gaussian, divergenceless field [37, 38] . 
For this reason, the set of ground states is said to constitute a Coulomb phase. 

The possibility that spin correlations, averaged over ground states, have a 
long-range component, is not self-evident. Indeed, one might expect the fact 
that there are a macroscopic number of ground state degrees of freedom to 
signal the opposite, since their existence implies that the set of ground states 
includes local degrees of freedom that can fluctuate independently. In turns 
out, however, that some ground state correlations are impervious to all local 
fluctuations: in this sense they can be said to be topological. 

The simplest way to appreciate the existence of long-range correlations 
within ground states is to start from p.2[) and the fact that the total spin L a 
of each frustrated cluster a vanishes within all ground states. A consequence 
of this on the pyrochlore lattice can be visualised with reference to Fig. 11.21 In 
particular, consider for any ground state the total magnetisation m(z) of the 
lowest plane of sites in this figure. Its value (which may be small - for example, 
of order the square root of the number of sites within the plane, if spins are 
randomly orientated within the layer) is perfectly correlated with the mag- 
netisation of other parallel planes making up the lattice. Indeed, let m(z + 1) 
be the magnetisation of the plane neighbouring the lowest one. Since the sites 
in both planes taken together make up a layer of complete tetrahedra, the 
overall magnetisation of the layer is zero, and so m(z + 1) = — m(z). By ex- 
tension, m(z+n) = (—1)™ x m(z) for any n, a signal of long range correlations. 
The correlations give rise to sharp features, termed pinch points or bow ties 
in the Fourier transform of the two-point correlation function, averaged over 
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ground states, as obtained from simulations [39], large- n calculations [40], and 
diffuse neutron scattering measurements [41]. These singularities distinguish 
the diffraction pattern of the frustrated system from that of a paramagnet, 
but are weaker than those of Bragg peaks arising from Neel order. While their 
structure can be understood by building on our discussion of m(z) [10], a more 
complete approach uses a long-wavelength description of ground states. 

This continuum description of the Coulomb phase is obtained by mapping 
spin configurations onto configurations of vector fields in such a way that the 
ground state condition L Q = 0, involving the specifics of the lattice structure, 
is translated into the requirement that the vector fields have lattice divergence 
zero. This second version of the constraint has the advantage that it can be 
implemented in the continuum [37,38]. To describe the mapping in detail, it 
is necessary first to discuss some features of the simplex lattice, introduced 
in section 11.21 We require the simplex lattice to be bipartite. This is the 
case, for example, for the diamond lattice, the simplex lattice associated with 
the pyrochlore lattice. (In models without a bipartite simplex lattice, the 
correlations discussed in this section are absent [42].) For a bipartite simplex 
lattice, one can adopt an orientation convention for links, taking them to be 
directed from the simplices of one sublattice to those of the other. In this way 
one can define unit vectors ej oriented according to the convention and placed 
at the mid-points of links of the simplex lattice, which are also the locations 
of spins. Considering in the first instance Ising spins, a spin configuration is 
represented as a vector field B(r) on the lattice via 



The condition satisfied in ground states, ^2 i£a Si — 0, fixes the lattice di- 
vergence of B(r) to be zero. More generally, for n-component spins Si we 
require n flavours of vector field B ; (r), with I — 1, 2 ... n. These are related 
to spin components S\ via B'(r) = S\ei, so that in ground states each flavour 
has divergence zero. Note that the fields B'(r) are defined in real space, and 
that the global 0(n) symmetry of the spin Hamiltonian is implemented as a 
transformation within the space of flavours, I. 

Continuum versions of these vector fields B'(r) result from coarse-graining 
and the restriction to ground state configurations is imposed exactly by requir- 
ing the continuum fields to be solenoidal. Each coarse-grained state represents 
many microscopic configurations and should have an entropic weight that re- 
flects this. It is plausible that small continuum field strengths will arise from 
many different microscopic configurations, and that large field strengths will 
be generated by fewer microscopic states. This suggests [37, 42] the weight 



This theory has a single parameter, the stiffness k, whose value affects the 
amplitude but not the form of correlations and is determined microscopically. 



B(rj) = Siii . 



(1.11) 




(1.12) 



1 Fluctuations in geometrically frustrated magnets 



15 



In d = 3 dimensions all other terms consistent with symmetry that might be 
added to the effective action are irrelevant in the scaling sense, and so (|1.12|) 
is expected to have a universal validity. The resulting correlation function, 

(BH0)Bp (r)> = |=L (■'"'■'> y 2 *") , (1.13) 

falls off with a fixed, integer power of distance, and has a characteristic, dipolar 
angle dependence. 

The fixed, integer power appearing in (|1.13|1 stands in contrast to be- 
haviour in two other situations in statistical mechanics for which power-law 
correlations appear: those of a system undergoing a continuous phase transi- 
tion, and of the low-temperature phase in the xy model. The form of correla- 
tions in (| 1 . 1 3(1 is instead similar to those generated by Goldstone modes in the 
ordered phase of a system with a spontaneously broken continuous symmetry, 
and an equivalence between that and the Coulomb phase can be developed 
by passing to a dual description of the frustrated magnet [43]. 

At finite temperature thermal fluctuations out of the ground state manifold 
generate a finite correlation length £ which acts as a cut-off for the power- 
law in (|1.13p . The scale £ diverges at low temperature, as £ ~ T -1 / 2 in a 
Heisenberg model, and exponentially in an Ising model. 

1.7 Dynamics 

As we have seen in some detail, geometrically frustrated magnets in the tem- 
perature window T c C T « lew are strongly correlated, yet lack long range 
order. Their dynamics in this regime has novel features which we summarise 
in this section. 

An obvious first step to understanding low temperature dynamics is to 
apply harmonic spinwave theory, starting from one of the ground states. The 
results one expects are summarised in terms of the density of states p(u>) in 
frequency u> in Fig. 11.101 Within the harmonic approximation, excitations are 
of two types. One type, similar to those in conventional magnets, forms a band 
of finite-frequency states, with a maximum frequency ~ 0(JS/h). The other 
type (those associated with the ground state coordinates x, in the discussion 
of (jl.lOp ) are zero modes. For example, for excitations from a generic ground 
state of the Heisenberg antiferromagnet on the pyrochlore lattice, one quarter 
are zero modes. 

There is a clear interest in understanding in more detail the nature of 
these zero modes. In the most cases, however there is a obstacle to a simple, 
analytical treatment, which stems from the fact that spin configurations in 
representative classical ground states do not have long range order. This means 
that, even though the lattice itself is periodic, the equations of motion cannot 
be diagonalised by Fourier transform, and results of the kind sketched in 
Fig. ll.lOl can be obtained only numerically. 
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p(co) 



zero modes 




spin 
waves 



Fig. 1.10. Density of states in frequency for harmonic excitations in a geometrically 
frustrated antiferromagnet 



To circumvent this difficulty and illustrate in a simple fashion how a disper- 
sionless band of modes can arise, it is interesting to consider a geometrically 
frustrated Heisenberg antiferromagnet in the presence of a magnetic field h 
strong enough that the ground state is fully polarised. Using the standard 
Holstein-Primakoff transformation to write operators for spin components S\ 
in terms of boson creation and annihilation operators a\ and a* , with 



Sf = S 



ST 



(25) 
(25) 



1/2 «, 



(1.14) 



we have 



a 3 +a 3 a i 



0{S°) 



(1.15) 

The right-hand form of (|1.15[) is a tight-binding model for bosons moving 
on the lattice with a nearest-neighbour hopping amplitude J5 and a chemical 
potential /i = zJS—h that is linear in the magnetic field h. It is a characteristic 
of the lattices we are concerned with that such a tight-binding model has a 
dispersionless band with eigenvalue — 2J5, which lies at the bottom of the 
spectrum for J > 0. 

Eigenvectors of the tight-binding Hamiltonian from the dispersionless band 
are straightforward to picture. One for the kagome lattice is represented in 
Fig. 11.111 Here, eigenvector amplitudes are zero at all sites except for those 
around one hexagon of the lattice, on which they have equal magnitude and 
alternating signs. The state is an eigenvector because there is destructive in- 
terference between hopping processes that move the boson off the occupied 
hexagon. It belongs to a dispersionless band since it is degenerate with many 
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other, equivalent states, based on the other hexagons of the lattice. The condi- 
tion for an arbitrary vector, with site amplitudes ipi, to be a linear superposi- 
tion of such states is that J^iea fa = I0r each triangle a. Both the extension 
of this condition to other lattices constructed from corner-sharing simplices 
and its parallel with the ground state condition in spin models, L Q = for all 
a, are obvious. There is a gap to excitations for large h, when fj, is large and 
negative. The gap falls to zero at the critical field strength h c — JS(z + 2) 
at which fi crosses the energy of the dispersionless magnon band. For h < h c 
these states are populated, and the magnetisation deviates from its saturated 
value. In this field range we recover the classical ground state degeneracy of 
the zero-field problem. 




Fig. 1.11. A magnon mode from the dispersionless band 

We now turn to a discussion of dynamics beyond the harmonic approxi- 
mation. As a starting point we should consider the full equation of motion 

1Q J 

h ^r = St x Hj = 2 {[La + Lp] x St ~ Si x [Lq + L/3]} • (L16) 

Here is the exchange field acting at site i. It is parallel to Si in a clas- 
sical ground state, but in excited states has transverse components. These 
can be expressed as shown, in terms of the total magnetisations L Q and 
of the clusters a and j3 that share the site i. For large S we can treat this 
equation of motion classically. Within the harmonic approximation, obtained 
by linearising the right-hand side, the exchange field is a superposition of 
contributions from finite frequency modes, which maintain phase coherence 
indefinitely, and therefore average to zero over long times. Anharmonic inter- 
actions have two consequences, which are distinct but turn out to be closely 
linked [9, 10]. One is to generate a lifetime for the finite frequency modes. The 
other is to introduce coupling between the ground state coordinates x and the 
coordinates y orthogonal to the ground state manifold. Over long timescales 
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this coupling drives the system around the ground state manifold. Within 
the framework of (| 1 . , this long-time component to the dynamics arises be- 
cause, once spinwaves have a finite lifetime, the exchange field is no longer a 
superposition of exactly harmonic contributions. Instead, on timescales longer 
than the lifetime, it is better thought of as a stochastic quantity. In turn, the 
long-time motion of the system around the ground state manifold is itself 
a source of dephasing for finite frequency excitations. Specifically, since the 
harmonic Hamiltonian is time-dependent, an adiabatic approximation is not 
exact, and modes are mixed at long times. There is a separation of timescales, 
since typical spinwave periods are fixed, while spinwave lifetimes diverge as 
T -1 / 2 and the timescale for motion between groundstates diverges faster, as 
T- 1 [10]. 

These ideas suggest a much simpler approach to calculating the spin auto- 
correlation function, in which we treat the exchange field as a stochastic quan- 
tity, using the equation of motion 

? = S(i)xH(t) (1.17) 

with the correlation function 

(H l (t)H m (t'))=rS lm 6(t-t') (1.18) 

for components H (t) of H(t). The noise intensity r can be estimated using 
equipartition. We have 

(Hi(0)-H((t)>dt ~ / (L a (0)-L a (i))dt. (1.19) 

-oo J — oo 

In addition, L Q (i) is a superposition of contributions with amplitudes from 
thermally excited spinwaves: 

L a (t) = J2^- iwit ■ (1-20) 
i 

Combining these assumptions, we find [10] 

r~ <|a w |>m|^ ~^. (i.2i) 

The Langevin equation (| 1 . 18|) itself is straightforward to solve, and yields 

(S(0) • S(t)) = S(S + 1) exp(-ck B Tt/hS) , (1.22) 

where c ~ 0(1) is an undetermined numerical constant. This result is notable 
for the fact that temperature alone sets the time scale: J drops out of the 
long-time dynamics in the low temperature regime. In this sense, behaviour 
matches that expected at a quantum critical point, although the underlying 
physics is quite different. 
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The predictions of Q1.22p have been tested both in simulations and in ex- 
periment. Molecular dynamics simulations proceed by direct integration of 
the equations of motion, (|1.16p . with an initial configuration drawn from a 
thermal distribution and generated via Monte Carlo simulation. Results for 
the classical pyrochlore Heisenberg antiferromagnet [9,10] in the temperature 
range keT <C JS 2 reproduce both the functional form of (| 1 .22(1 for the time 
dependence of the autocorrelation function and the scaling of relaxation rate 
with temperature. In experiment, inelastic neutron scattering offers direct ac- 
cess to spin dynamics. Early measurements of the energy width of quasiclastic 
scattering in CsNiFeg [44] were fitted to a Lorentzian, the transform of the 
time-dependence given in (|1.22p . yielding a relaxation rate for T < |@cw| 
that is strongly temperature dependent, although without specific evidence 
for the (subsequently proposed) linear variation with T. More detailed data 
for SCGO [45] confirms a relaxation rate of order ksT/h, and very recent 
measurements on Y2RU2O7 [46] display rather clearly a relaxation rate pro- 
portional to temperature. For the future, it would be interesting to have both 
theoretical and experimental information not simply on the autocorrelation 
function, but on dynamics as a function of wavevector. 

1.8 Final Remarks 

In conclusion, we have seen how geometrical frustration in classical magnets 
can lead to macroscopic ground state degeneracy and the suppression of long 
range order. Low temperature states in model systems, although disordered, 
are very different from those of a non-interacting paramagnet: correlations 
are power-law in space, and decay in time at a rate set by temperature alone. 
Many experimental systems display these features within the temperature 
window T c <T < |6>cw| where behaviour is dominated by nearest neighbour 
exchange. Behaviour in this regime is well summed up in the term coined by 
Jacques Villain [8], a pioneer in the field: Cooperative Paramagnetism. 
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